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Or: "How to generate an ensemble in a single event?" Following re- 



^ ■ cent work on entropy in strong interactions, I explain the concept of 

environment-induced quantum decoherence in elementary quantum me- 

• ^H , 

r> I chanics. The classically chaotic inverted oscillator becomes partially de- 

coherent already in the environment of a single other oscillator perform- 
ing only vacuum fluctuations. One finds exponential entropy growth 
in the subsystem with a Lyapunov exponent, which approaches the 
classical one for weak coupling. 
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Recently the long-standing "entropy puzzle" of high-multiplicity events in strong 
interactions at ultra-relativistic energies has been analysed from a new point of view 
HI]. This is related to the concepts of an open quantum system and environment- 
induced quantum decoherence. The problem dates back to Fermi and Landau and is 
intimately connected to understanding the rapid thermalization of high energy density 
(^ 1 GeV/fm ) matter P|. Why do thermal models work so well in reproducing 
global features of hadronic multiparticle final states? Why do they work at all? 

Or, Why does high-energy scattering of pure initial states lend itself to a statisti- 
cal description characterized by a large apparent entropy from a mixed-state density 
matrix describing intermediate stages in a space-time picture of parton evolution? 
Effectively, unitary time evolution of the observable part of the system breaks down 
in the transition from a quantum mechanically pure initial state to a highly impure 
(more or less thermal) high-multiplicity final state. Note that the unitary time evolu- 
tion operator, exp(— iHt), always transforms a pure state into a pure state, according 
to the Schrodinger equation, which cannot produce entropy under any circumstances 
(cf . below) . This was discussed in detail in Refs. [Q , where more references concern- 
ing formal aspects of this work can be found. Based on analogies with studies of the 
quantum measurement process ( "collapse of the wave function" ) [|| and motivated 
by related problems in quantum cosmology and by non-unitary non-equilibrium evo- 
lution resulting in string theory |^, I argued that environment-induced quantum 
decoherence solves the entropy puzzle of strong interactions. 

A complex pure-state quantum system can show a quasi-classical behaviour, i.e. 
an impure density (sub)matrix together with decoherence of the associated pointer 
states in an observable subsystem |]l], |^, Q. I will demonstrate in the following that 
the decoherence process is uniquely correlated with entropy production. Considering 
strong interactions, in particular, there is a natural Momentum Space Mode Separa- 
tion due to confinement, which is defined in the frame of initial conditions for the 
time evolution and for the physical (gauge) field degrees of freedom. Thus, almost 



constant QCD field configurations form an unobservable environment, which inter- 
acts with the observable subsystem composed of partons. The environment modes 
are unobservable, since they can neither hadronize nor initiate hard scattering among 
themselves, whereas the partons are observable in the sense of parton-hadron duality 
or deep-inelastic scattering; equivalently, low-energy coloured vacuum fluctuations 
cannot propagate into asymptotic states. 

Previously, I studied the induced quantum decoherence and entropy production 
in a non-relativistic single-particle model resembling an electron coupled to the quan- 
tized electromagnetic field, however, with a deliberately enhanced oscillator spectral 
density in the infrared. The Feynman- Vernon influence functional technique for quan- 
tum Brownian motion provided the remarkable result that in the short-time strong- 
coupling limit the model parton behaves like a classical particle |l|]: Gaussian parton 
wave packets experience friction and localization, i.e. no quantum mechanical spread- 
ing, and their coherent superpositions decohere. The decoherence process has been 
shown to lead to entropy production in this oversimplified parton model. 

It seems somewhat more realistic to consider two coupled scalar fields represent- 
ing partons and their non-perturbative environment, respectively. In the functional 
Schrodinger picture employing Dirac's time-dependent variational principle, i.e. a 
non-perturbative method, I derived a Cornwall- Jackiw-Tomboulis (CJT) type effec- 
tive action and the equations of motion for renormalizable interactions |]T|. Thus, 
analysis of the entropy puzzle in strong interactions leads to study an observable field 
(open subsystem) interacting with a dynamically hidden one (unobservable environ- 
ment), i.e. quantum field Brownian motion. 

Summarizing, my point of view is that partons feel an unobservable (gluonic) 

environment, which manifests its strong non-perturbative interactions on a short time 

scale (-C 1 fm/c) through decoherence of suitable partonic pointer states 0, their 

■^In general, these are not single-particle states but rather coherent (Gaussian) wave functionals, 

as constructed in the second of Refs. ||l|. 



quasi-classical behaviour, and entropy production. If confirmed in QCD, this will have 
important consequences for parton-model applications to complex hadronic or nuclear 
reactions. The emergence of structure functions from initial-state wave functions can 
and will be further studied in this approach. 

Instead of representing the formalism and more technical results from Refs. [Q, 
I want to demonstrate here in simple quantum mechanical examples the basic Why 
and How of the solution to the entropy puzzle. 

Consider a system that can be described in terms of two normalized discrete 
basisstates, |1) and |2). Forming a pure state, \ip) = ai|l) + a2|2), by a coherent 
superposition with amplitudes ai = p^^"^ and 02 = {I — pY^"^, the corresponding 
density matrix, p = \iIj){iIj\, is 
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where p^ is obtained by diagonalization. Note the off-diagonal interference terms in 
Pij. Furthermore, observe that p^ has only one non- vanishing eigenvalue. Introducing 
the von Neumann or statistical entropy, 

S[p] ^ -Tiplnp , (2) 

we find S[p] = S[p^] = —[1 In 1 + In 0] = 0, i.e. no entropy in a pure state. Secondly, 
forming a mixed state (ensemble) such that the system is in state |1) with probability 
p and in state |2) with probability 1 — p, the density matrix becomes p' = |l)p(l| + 
|2)(1 — p){2\, i.e. a decoherent superposition. Hence, we obtain 

( p \ 

P[, = • (3) 

\0 l-pj 

The density matrix (^ shows no interference terms and is diagonal per se.0 Then, 
S[p'] = S{p) = —[p\np+ (1 — p) ln(l — p)] 7^ 0, generally. In fact, < S{p) < 
""^Note that Tr p = Tr p' = 1: the system is in some state with total probabihty 1. 



S{l/2) = In 2. Total ignorance about the state of the system {p = 1/2) corresponds 
to In 2 units of entropy in a two-state system, i.e. 1 bit of information is lost compared 
to certainty about its state {p = 0,1). 

One concludes that entropy production can only occur if the interference terms 
of the density matrix representing a more or less pure state of the observed system 
decay dynamically.^ In a closed system evolving unitarily in time, however, there is 
no way to transform, for example, p into p', see Eqs. (||), (^. Only the interaction 
of the system with an environment |T|, ^, ^, can have such an effect. The presence 
(and integrating out) of the environment degrees of freedom essentially changes the 
dynamics of the observed system. This can lead to the decay of the interference 
terms in its density matrix, i.e. environment-induced quantum decoherence, which is 
necessary to increase its impurity and, thus, to produce entropy. 

Next, consider a non-relativistic particle moving in a one-dimensional double- 
well potential presenting the observable subsystem, which is coupled translationally 
invariant to a single environment oscillator. The classical action is 

S = I dt [\Mx^ + \my^ - \muj^{y - xf + \M^^x^ - ^yAMfiV} . (4) 

For simplicity, lei M = m = VL. Then, properly rescaling by f2, one obtains 

S = j dt [\x^ + \y^ - \uj\y - xf + ix^ - lAx^} , (5) 

in terms of dimensionless quantities and two coupling constants, cu^, A > 0. For cu = 

the minima of the doublewell lie at x± = ±(A/3!)^^/^, at a depth of — 3/2A (the local 

maximum is zero at a; = 0). Presently, I want to study the case that the excitation 

energy of the x-particle (X) is smaller than the level spacing uj of the environment 

oscillator (Y), which is assumed to be in its ground state. Starting with a given initial 

state of X, I will calculate the time evolution of the corresponding density matrix px 
^The argument does not depend on particular physical characteristics of the system; it holds for 
the two-state system as well as for an interacting quantum field. 



under the influence of tlie vacuum fluctuations of Y. Excited states of Y contribute 
only virtually here; they cannot become real due to energy conservation. 

For illustration, I choose the metastable initial state of X, when classically the par- 
ticle "rests on top of the hill" {x = 0). Quantum mechanically this can be represented 
by a minimum uncertainty Gaussian wave packet, 

^{x,t = 0) = 7r-i/V'^Vi^V< , (6) 

with wq <^ A^^/^. Also, assume to to be sufficiently larger than 3/2A. 

First of all, let the system evolve classically. Nothing will happen. However, any 
infinitesimal perturbation of the fine-tuned initial conditions causes X to move "down 
the hill", left or right (L or R), dragging Y along. There is local chaos in the sense of 
extreme sensitivity to the initial conditions at x, y ~ 0; arbitrarily small uncertainties 
in the initial conditions lead to a loss of predictability. For an ensemble of initial 
conditions X switches with probabilities piit') and PB.{t') = 1 — PL{t') between L and 
R, respectively, if at least one trajectory passes x = in a certain interval [f — e, t' + e]. 
This corresponds to a loss of information about the actual binary decision "either L 
or R" and an entropy Sx = —J2i=L.RPi^^Pi-f\ Note that Sx or pi^n are strongly 
conditional ("fine-grained") quantities. In distinction, the usual classical entropy is 
calculated after "coarse graining", i.e. by constructing a local probability density 
f{t) in the phase space of X related to the ensemble average over initial conditions, 
Sc.g.{t) = J dxdp f(t) In fit) + C. A chaotic loss of predictability from strongly 
diverging trajectories in phase space causes Sc.g. to increase: as time passes, more 
and more cells of the coarse graining contribute — effectively, the phase space volume 
occupied by the ensemble grows (without violating Liouville's theorem). 

In conclusion, in a classical system, be it chaotic or not (with or without coarse 

graining), entropy can only be produced IF there is a physically relevant ensemble of 
^If the ensemble of initial conditions is constrained to preserve the reflection symmetry of the 
action, Eq. (|[), then PL{t') —PR{t') — \ and Sx{t') — In 2 stay constant. 



initial conditions. Thus, one cannot explain altogether classically entropy production 
or thermalization in a single high-multiplicity event in strong interactions. 

Secondly, let the system evolve quantum mechanically. To begin with, let 
there be no coupling to the environment {u = 0). Even with the fine-tuned initial 
condition, Eq. (^), the amplitude ip to find X at a particular space-time point begins 
to fiow "down the hill" symmetrically (L and R) due to the quantum spreading of 
the wave packet. For a free particle wq — ^ w{t) = (wq + W7(7^t^)^/^ (for M = 1); here 
one expects an accelerated spreading "downhill", cf. Eq. (jT^) below. The related 
probability density \ip\'^ also evolves and stays symmetric; generally, it cannot be 
simulated by the classical evolution starting with an ensemble of initial conditions due 
to the absence of quantum interference between classical trajectories. In any case, the 
system remains in a pure quantum state. The density matrix is px{t) = \ip{t)){tlj{t)\, 
{ipit)) = exp[— 'iHot]|'?/'(0)), where Hq is the Hamiltonian of X from Eq. (P) with uj = 0. 
Therefore, S[px{t)] = — llnl = 0, cf. Eq. (^. Quantum mechanically one knows 
ever5d;hing there is to know about a closed system (X), given any pure initial state 
and its Hamiltonian, which consistently yields 5'[px] = 0. Even with an ensemble of 
initial states, i.e. an impure density matrix px(0), there is no entropy production, 
since S[px(t)] = S'[exp(— 'iHo^)px(O) exp(+zHo^)] = S[px{0)] stays constant.[] 

The situation changes completely if the subsystem (X) evolves quantum mechan- 
ically coupled to the vacuum fluctuations and virtual excitations of the environment 
(Y). With the above assumptions the initial density matrix of the total system is: 

p(t = 0) = Px{0)®Py{0) , (7) 

with matrix elements px (a;, x'; 0) = tv^^^'^Wq^ exp[—^{x'^ + x''^)/wq] and py(y, y'; 0) = 

(cij/vr)^'^ exp[— |a;(t/^+l/'^)]. The time evolution of the density matrix of the observable 

subsystem, pxit) = Tryp(t), can be calculated with the Feynman- Vernon infiuence 

functional technique; I will make use of general results obtained in the first of Refs. 
^The unitary (time evolution) transformation docs not change the eigenvalues of px- Thus, the 
statistical entropy, Eq. (^ , cannot possibly show a sign of classical chaos in a closed system. 
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m. The idea is to derive a propagator for px-, which incorporates the influence of the 
environment degrees of freedom (Y) exactly. This can be achieved, since Y and its 
couphng to X are at most quadratic in coordinates and momenta, see Eq. (^. 

It should be remarked that the final state of the environment is not specifled; 
presently, it may contain virtual excitations of Y.^ The relevance of this for a high- 
multiplicity hadronic (or nuclear) reaction is the following: Even though the QCD 
vacuum "far away" conforms to the usual one before and after, the additionally 
produced secondary hadrons all require a dressing of their valence quarks by localized 
virtual excitations of the vacuum or environment, which obviously makes an essential 
difference as compared to the initial state. 

Presently, the resulting density matrix px{t) is (cf. also the flrst of Refs. [|I[]): 

px{z-,z^,t) = 7r-i/2«;-i(t)e-[^+-^W^]V^'(i) 

X e-^-{^ + iWc' - d-'[B + lwibc]yw\t)} 

^ ^iz.{{az+ - 2d-'[B + K6c][z+ - t;(t)t]/«;2(t)} ^ ^g^ 

with the effective velocity v{t) = for the zero-momentum initial wave packet, the 
effective width w{t) = 2^\d\~^, ^ = {A + ^Wq"^ + jW^b'^Y^'^, and with rather com- 
plicated time-dependent coefficients A, B, C, a, b, c, d, to be discussed elsewhere; 
the coordinates in Eq. (|^) are z_ = x ~ x' and z+ = ^(x + x') in terms of ordinary 
one-dimensional ones. Since we are particularly interested in the decoherence process 
and entropy production, we consider only the simplest off- diagonal density matrix 
elements here, px{x,x' = —x,t) = px{z- = 2x,z+ = 0,t). They can be directly 
related to the linear entropy produced in the observable subsystem (X): 

/OO /"OO 

dz. dz+pxiz_,z+,t)px{-z.,z+,t) 
-OO J —OO 

= 1 - Ww-^ , (9) 



®As a corollary to the Schmidt decomposition 11] it is easy to prove that starting with an overall 
pure state of the complex system, cf. Eq. (M), IF the final state of the environment is a pure state, 
THEN the observable subsystem ends up in a pure state too (without entropy production). 



with ci = C + jWqC^ — {B + ^WQbc)'^/{dwy and independently of the initial wave 
packet momentum (p = at present). Thus, inserting (^ into (|]), one obtains: 

px{x,-x,t) = 7T-^/^w~\t) exp{-x^w-\t)[l-S''''{t)]-^} , (10) 

dxpx{x,-x,t) = l-^'^"(t) > e-^W . (11) 



The inequality results from the fact that the linear entropy provides a lower bound for 
the relevant statistical entropy, cf. Eq. (^, as shown in |jl|. Note that Eqs. (P) — (|lTl) 
are completely independent of the time- dependent functions entering there, which are 
specific for a particular dynamical system. They are based, however, on the Gaussian 
structure of the subsystem density matrix, Eq. (^. 

At this point the attentive reader might wonder what happened to the non-linear 
interaction oc Ax^ of the double- well potential, see Eq. (|5|). Of course, it cannot 
be treated exactly. I employed a mean-field- type approximation, replacing ^Ax^ by 
^X{x^)x^ = |A^(t)x^. As long as one studies only the initial time-evolution over short 
periods, as compared to the time a classical particle would need to "roll down the hill" , 
one may even set A ^ 0. For the following qualitative considerations, A plays the role 
of an adiabatically changing parameter. However, a more accurate approximation is 
necessary (and feasible) to follow the truly long-time quasi-periodic motions of the 
system. Then, one expects periods of increasing decoherence and entropy production, 
cf. below, followed by periods of quantum revival in the observed subsystem. The 
more complex the environment becomes, the more unlikely quantum revivals will be, 
since the total system including the environment finds more and more ways to evolve 
before a reconstruction of the subsystem initial-state wave function.^ 

To begin with, it can be checked explicitly that there is no entropy production for 

a vanishing coupling to the environment, S'^*!!o(t) =0, cf. Eq. (|]). Next, calculating 
^Such effects have been experimentally observed in even simpler systems involving a two-state 
subsystem of one Rydberg atom coupled to a single mode of the electromagnetic field . 



the effective width in the long-time limit [^, one finds: 

w{t) = (W + ^o'f=' + -^'f=V' + u^r^y^' yr^ e^P t- > (12) 

witht_ = f_t, /± = [±icJ+2^(icu^-cuV^)l/2]V^ andu;J = ±uj^ + A'^~l. Assuming 
a sufficiently small coupling, lu"^ < 1, note that tuj is negative as long as A^(t) < 
1 — uj'^. Thus, the width grows exponentially with an effective Lyapunov exponent 
/_. For vanishing coupling to the environment, it reduces to the classical Lyapunov 
exponent of the inverted oscillator, /^^° = (1 — A^)^/^, while the width becomes 
Wu,=o{t) = {wq + [wo/-"°]"^)"'^^^exp(/^=°t). This suggests quite generally that the 
time-dependent widths of suitable (Gaussian) wave packets may serve as "quantum 
indicators" of chaotic behaviour in the corresponding classical system. 

It is remarkable how the Lyapunov exponent reffects the dynamics: as the wave 
packet spreads "downhill", A^(t) oc (x^) increases until /_ reaches zero (becoming 
purely imaginary afterwards), when uj_ = 0. At this point the behaviour becomes 
regular in the sense of being governed by harmonic motions close to the minima of 
the double-well potential with a correspondingly milder time- dependence of the width 
(cf. the model studied in the first of Refs. ^). 

The second dynamical time scale f^^ always stays real. It is relevant for certain 
non-Markovian effects generated by the interaction with the environment (/^^° = 0). 
These become clearly visible in the entropy evaluated in the same limit as Eq. (pT^). 
Using ci = C + 0{wq), one obtains in leading order: 

with t+ = /+t. Thus, the linear entropy approaches exponentially its saturation 

value 1 on the time scale set by the Lyapunov exponent, see Eq. (|12[), and the 
^"^Here, hyperbolic functions dominate over trigonometric ones in the time-dependent coefficients 
in Eq. M) for times such that a classical particle would still be "rolling down the hill" of the 
potential; this restriction is presently assumed for simplicity. 







(13) 


2u;^cost+ 2 2r^^'^^+ 


COSt+,2 
/- J 


-1/2 



von Neumann entropy grows exponentially according to eq. (pT]), at least as fast. 
Note that the periodic function multiplying w~^{t) in Eq. (^) is approximately 
oc cu-i/2[l + 2^/2 + sin2(2i/4cut)]-V2 for /_ ^ /+ and oc 07-^2 [l + cos\2'/''ujt)]-^/^ 
for /_ ^ 0; i.e. it persists qualitatively even until the effective Lyapunov exponent 
becomes imaginary, when the stabilizing effect of the x^-term in the potential is felt. 
To summarize, in the above specified long-time limit and for the chosen initial 
conditions, Eq. (0), one obtains the observable subsystem (X) density matrix, 

px{z-,z-,,t) = vr-/2^-(t)e-W + i^-[l-^'"W]"'}^"'W + --+/- . (14) 

Even though this density matrix describes the exponential entropy production and 
its diagonal matrix elements with z_ = grow rapidly, apart from the overall nor- 
malization factor, the (simplest) off-diagonal matrix elements {z+ = 0) do not really 
decay here as in usual models of quantum decoherence |[l|, Q. This is no surprise in 
view of the "poor environment" considered at present, which has only one degree of 
freedom frozen in its ground state (modulo virtual excitations). Loosely speaking, 
it is unable to accommodate all the phase information contained in the off-diagonal 
density matrix elements of the subsystem. 

In conclusion, a strong observable entropy production in a quantum system, which 
shows a chaotic behaviour in the classical limit with exponentially growing modes, re- 
quires only a minimal decohering effect due to an environment of vacuum fluctuations 
coupled to it from a higher energy scale. 

In particular, complementary to previous studies 0, P], one observes here that 
the environment does not necessarily have to be at any finite temperature for this 
effect of partial decoherence to work. Furthermore, it should be realized that the 
Schmidt decomposition reveals the remarkable fact that the density submatrices of 
"subsystem" and "environment" always have identical non-zero eigenvalues [|l| . Thus, 
from the point of view of calculating the entropy, Eq. (^, their roles can be inter- 
changed and it is a matter of practicability to decide which part of the total Hilbert 
space is integrated out to find the entropy of the physically observed subsystem. The 

10 



environment-induced quantum decoherence and its relevance for entropy production 
have presenly been illustrated by an elementary example, which, however, points out 
to interesting consequences for the quantum evolution of classically chaotic non-linear 
field theories. 

I thank P. Carruthers, M. Danos, N. E. Mavromatos, B. Miiller, O. Bertolami, 
and J. Rafelski for stimulating criticism and helpful discussions. 
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